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We report on the experimental observation of weak localization in an optically 
induced disordered (2+l)-dimensional photonic structure. Our flexible method of 
optical induction is applied with a nondiffracting random intensity distribution. We 
focus on the analysis of a statistical output spectrum for many probe events with 
variance of the incoming beam’s transverse spatial frequency. For particular spatial 
frequencies we find considerable signatures of transverse coherent backscattering. 


Coherent backscattering (CBS) is a pure wave effect evolving from scattering in random 
media. Its appearance can briefly be described as incoming wave vectors are in average reversed 
such that the probability to find outcoming waves in the antiparallel direction of incidence is 
enhanced. |1| This effect bases upon wave interference and, given that a random medium has 
appropriate conditions, it can be found among various systems where sufficient long wave 
propagation is provided. Thus, numerous experiments were designed and realized in order 
to explore the conditions yielding to CBS, and to prove the universality of this effect, e.g. 
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as well as matter |6,[7] and light waves in turbid 
Aiming to find a wave which scatters preferably 


by applying acoustic |2-4| and seismic 
colloids |8j[9] as well as in cold atoms JTc 
elastically, especially by using light as the preferred wave, a lot of efforts have been invested 
during the last three decades in order to clarify scattering processes inside a turbid medium 
12j[l3—a comprehensive review can be found in Ref. 
which is in consistence with all experimental results 
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and to establish a universal theory 
16 . Since the scaling theory predicts 


that there is a transition from extended to localized states in a system of more than two 
dimensions, 17 special attention has been payed to describe localization in three-dimensional 
samples. 14] However, from experiments on three-dimensional random media it is still hard 


to distinguish whether reduced transmission or reflection rates come from localization or from 


absorbing and resonant processes. {18,19 Still, especially identifying and examining the critical 
regime between localizing and diffusive behavior is a burning issue to be answered. 1201 


Though all modes are localized in one- and two-dimensional random systems, media of these 
dimensionalities are of high interest to investigate localization and backscattering events in 
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various configurations as they allow for a deep insight into these wave events in a less complex 
manner. 0 21 For there is a distinction between strong and weak localization (commonly 
referred to Anderson localization 1221 and Coherent Backscattering), it is moreover most im¬ 
portant to find model systems in which can be switched between both regimes and additionally 
compared to ballistic propagation at will. It has been shown in several publications 123 -26 


that in an optically induced randomized potential ensemble the observation of transverse strong 
localization is achievable by defining proper conditions. 

In this contribution we describe the preparation of proper ensembles of two-dimensional ran¬ 
dom photonic potential landscapes to experimentally observe and investigate CBS signatures. 
The potentials in terms of disordered refractive index modulations are optically induced in a 
photorefractive strontium barium niobate crystal (SBN). 

In order to induce such a configuration we employ so-called nondiffracting random inten¬ 
sity patterns (NRP) 128] with the transverse intensity being randomly distributed. To get an 
impression of these intensity distributions, a three-dimensional volume is scetched in Fig. [I] 
(a) where the trihedron clarifies the orientation and direction of propagation which is congru¬ 
ent with z. As for all nondiffracting beams, this distribution is constant along a particular 
propagation distance, i.e. the intensity is independent of the third dimension, according to 
the direction of propagation. 129] Selecting appropriate parameters for the structural size of 
such intensity—for these random beams we introduced the nomenclature photonic grain size 
c7 in an earlier publication 26 —one can reach nondiffracting distances of several centimeters 


and even longer. Throughout all the experiments presented here, we fix a to 20 pm which is 
connected to a nondiffracting distance that is longer than the long site of the SBN crystal of 
dimension 5 mm X 20 nun x 5 mm. The nondiffracting property makes these light configura¬ 
tions particularly appropriate to implement (2+1) dimensional systems and thus to mimic the 
temporal development in a two-dimensional photonic configuration. |30 

With these ingredients to form (2+l)-dimensional, transversely disordered potentials we 
establish an adequate platform to examine weak transverse localization, involving coherent 
backscattering events in transverse direction, i.e. perpendicular to the direction of light propa¬ 
gation. It has already been shown that regimes of strong localization can be created in order to 
present transverse Anderson localization in disordered regular but also completely randomized 
photonic lattice systems. 123,24,26| 

For the induction process of one single potential configuration we supply an external field of 
a few kV cm” 1 and a writing light field whose intensity contains the structure of the desired 
potential and which is ordinarily polarized with respect to the crystal’s symmetry axis (c-axis). 
During the induction process charge carriers of the SBN crystal rearrange according to the 
intensity modulation and cause internal electrical fields which in turn induce a refractive index 
modulation due to the linear electro-optic effect 31 . Due to the (2+l)-dimensional character 


of the potential, the induced refractive index distribution can be taken as an arrangement 
of waveguides where randomization is introduced with respect to the waveguide depth and 
position. Accordingly, every waveguide is potentially a scattering center with the ability to 
change the wave vector k of an incoming light field. 

The character of CBS can be analyzed best when averaging over a large ensemble of random 
potentials. Introducing an appropriate ensemble, we experimentally use 200 various NRP 
configurations, for the induction of each potential we apply an external field of 2kV cm -1 and 
set the illumination time of an individual writing beam to 10 s. In Fig. [I] (b), the numerically 
computed power spectrum of an effective intensity resulting from the average of 100 transverse 
NRP distributions is presented. Circular spectra of all contributing NRP have the spectral 
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radius k m , whereas the limiting spatial frequency for the power spectrum is 2k m . This can 
be understood by simple mathematic arguments as the power spectrum is the convolution of 
the NRP spectrum of radius kt with itself. As the spectrum of a NRP is a circle, the power 
spectrum of the effective intensity accordingly is a filled disk of radius 2 kt- For comparison to 
the mean intensity spectrum, a numerical calculation of the power spectrum resulting from the 
average of all random potentials is depicted in Fig. [l] (c). Due to an orientation anisotropy in 
SBN crystals, 32 the power spectrum of this mean refractive index modulation gets deformed 


as the modulation perpendicular to the c-axis drops to zero. This characteristic is prominent 
in Fig. [l] (c) where many scattering centers can be found in k x direction in contrast to k y 
showing no modulation just along this direction. By reasons of this anisotropy, we limit our 
experiments to the investigation of scattering characteristics in direction of k x . 

The probing process of each potential is performed in the linear propagation regime (cw 
laser light of a feq pWs at a wavelength of A = 532 nrn) with a broad Gaussian probe beam 
(GPB), polarized extraordinarily to the c-axis. Holding a fixed Gaussian beam waist of wq = 
200 pm that covers several scattering centers at the input facet, we introduce a variable tilt 
of the probe beam by a nonzero transverse component /q n in order to address various spatial 
input frequencies. All beams to optically induce the potential as well as to probe the induced 
potential are generated by use of a set of spatial light modulators 26 . Single output spectra 


were recorded by a CCD camera placed in the Fourier space of the SBN’s output facet. We 
further subsequently average all output spectra of induced potential realizations receiving a 
spectrum of probability P(ki n ) to identify most probable spatial output frequencies. To do 
so experimentally, we benefit from the reversibility property of the SBN crystal since after 
probing one potential configuration it can be erased by white-light LED arrays in order to 
induce another structure afterwards. After typical erasure times of 30 to 40 s, the refractive 
index distribution is homogenized again. 

In Fig. [2] a schematic illustrates the projection to the x-z plane of one of the considered 
(2+l)-d potential realizations where scattering centers are represented by gray straight bars. 
Four possible paths of light are indicated propagating from the input to the output facet of the 
potential, each with the same input wave vector k- m and individual output vectors k (mt . Path 
B thereby is the time inverted path of A, thus, A and B have the same lengths. Additionally, 
k[ n = — k ou t is valid (which is also valid for path D , but the path length is different from dor F). 
Hence, light fields that are propagating along A and B or along according path pairs interfere 
always constructively in the far field, since their transverse output frequency component is 
the negative of the input component. In contrast to this, light that propagates along other 
paths may interfere constructively, destructively or, in general, can have any modulus of phase 
difference between 0 and 27 t. In consequence, these interferences average out after considering 
many potential realizations. Solely, path pairs such as A and B will contribute constructively 
to the averaged output spectrum, namely exactly at the spatial frequency k ou t = — h m - 

The left column in Fig. [3]shows the input spectra of GPBs of various tilts from ki n /kt = 0.8 to 
1.2. For comparison, gray dashed circles illustrate the input spectrum of the potential-writing 
NRPs with fixed transverse wave vector component of length kt = 7r/o\ Corresponding prob¬ 
ability spectra after probing all realizations of the disordered potential are shown in the right 
column. Plots appended in the lower part of each image represent the azimuthal spectral 
distribution P{ip) at radius A’; n , where the angle ip equals zero at 12 o’clock and increases coun¬ 
terclockwise (cf. upper right image in Fig. Dj). In the right-column plots, input frequency k\ n 
(dashed green line) and expected CBS frequency fccBS (dashed red line) are marked. Following 
the spectra in the right column from top to bottom, we find that CBS signatures are present 
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for all tilts of GPB as the intensity is increased around the spatial frequencies &cbs where 
the CBS peak is expected, respectively. However, especially for a GPB tilt of k- m = 0.9 k t , we 
clearly see indications of CBS since for this parameter the /q n component of the mean output 
spectrum has the most significant contribution around the spectral position kcBS- 

To analyze this behavior more quantitatively, relative powers of particular spatial frequencies 
P{k) are plotted against varied tilts in Fig. [4j Here P(k ) is normalized to the overall power of 
the respective spectrum. For small tilts around k- m /kt = 0.8 the CBS peak P(kcBs) (blue data 
points) is higher than for the input frequency P(ki n ) (black data points). This behavior changes 
with increasing tilt as for high values of k\ n /kt there is almost no light backscattered to the 
specific expected CBS frequency. It is obvious that the scattering probability is weakened for 
higher spatial probing frequencies, as a considerable percentage of light is transmitted rather 
than scattered (cf. Fig. [3] spectrum for k\ n /k t = 1.2). 

Thus, since scattering processes are less probable for higher spatial probing frequencies, it 
seems that the propagation distance of 20 mm is not sufficiently long for a prominent CBS 
peak to be established. For these short distances the probability for light to find proper paths 
are rather small, however increase with longer distances. Consequently, by changing the tilt of 
our probing beam we can adjust the scattering strength of our potential. 

Worth mentioning is the texture of the probability spectra apart from fc; n and kcBS- Finding 
an almost filled spectrum for small tilts, the spectrum for higher k\ n values becomes more and 
more ring-like. The radius of this ring equals the spatial input frequency and spectral parts for 
small spatial frequencies become rather weak. Hence, for higher tilts the transverse wave vector 
component is almost conserved but merely changed in orientation in contrast to flat-angled 
probing below k in . For these cases, a higher percentage of light is scattered to the spectral 
center. Besides the direction, also the length of the transverse wave vector is changed. 

Elastic scattering processes can be described by applying the Ewald construction to our 
mean power spectrum of the potential that holds a limit frequency of 2kt (cf. Fig. [I]). Thereby 
all possible elastic scattering centers can be found on a circle of radius k\ n (the so called Ewald 
circle in two dimensions), touching the power spectral center which, together with the center 
of the Ewald circle, resembles h m . The amplitude along the latter circle thereby determines 
the scattering strength. Now it becomes obvious why scattering is weaker for higher k m —the 
Ewald circle includes decresingly strong scattering centers as k\ n increases (cf. Fig. 0 (c)). 
Thus at 2kt our potentials reveal an effective mobility edge, where, approaching this edge, 
scattering becomes more and more improbable. 


33,34 


However, coherent processes such as CBS cannot be explained with this model since the 
influence on the phase is neglected. In this connection, it is also not allegeable how processes 
that lead to a change of the transverse wave number k m (not only the direction) happen. Most 
likely, these processes go back to the influence of waveguiding for it is more probable that light 
is guided for a nonnegligible distance along local refractive index maxima before coupled to 
neighboring waveguides. This scheme is supposed to change the transverse spatial frequency 
markedly. 

In conclusion, we presented an optically induced disordered (2+l)-dimensional photonic 
system as a platform for the investigation of transverse coherent backscattering effects. The 
presented results reveal that our approach is explicitly appropriate to introduce potential 
ensembles of variable disorder strength where both regimes, strong and weak scattering, can 
be examined. We demonstrated that besides Anderson localization also CBS is observable. 
However, a few open points to be discussed are still left. The photonic grain size a and 
the angle of incidence k\ n turn out to be appropriate as control parameters for the strength 
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of localization. Though we previously analyzed the strength of disorder on the localization 
behavior 26], as well, the absolute value still needs to be identified and associated with the 
factor kl* |35|, including the mean free path length l* assumed to be in the range of a. This can 
be done by analysis of the CBS cone provided that a good signal to noise ratio is given. 36 


However, it is not clear so far if the shape of the cone in our case is fully consistent with 
the suggested models. Additionally, we identified a remarkable influence on the mean output 
spectrum by the shape of the refractive index power spectrum with upper limiting frequency 
k t representing an artificial, or effective, mobility edge. Can its existence possibly mimic the 
transition in three dimensions? 

All in all, subsequent experiments and analytics are necessary to characterize the presented 
system and to further prove its excellence for getting deeper insights into weak and strong 
localization phenomena. In particular, it would be interesting to find a way of breaking the 
time reversal symmetry by dephasing pulses to quantitatively identify the influence of coherence 
in our observations. [37] Another phenomenon whose experimental implementation is highly on 


demand is coherent forward scattering of light [38 39 where our approach could be a candidate 
to provide the proper conditions. 
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Figure 1. (a) Exemplifies a three-dimensional plot of volume of intensity for a single nondiffracting 
random intensity distribution. Modulation of intensity only occurs in x-y plane and agrees in every plane 
along direction of propagation z. (b) Power spectrum of effective intensity resulting from incoherent 
superposition of 100 different NRPs of spectral radius kt (position marked by dashed circle), (c) Power 
spectrum of corresponding averaged refractive index distribution originating from effective intensity. 
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Figure 2. Schematic of waveguides of the described (2+l)-dimensional photonic system of disorder. 
Gray lines represent projection of waveguides to the x-z plane. Lines are possible light paths scattering 
through waveguide arrangement. Blue path B is reversion of A (red), only those path pairs contribute 
to coherent backscattering. Gray paths symbolize arbitrary light paths with random phase difference 
to A or B. 
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Figure 3. Upper row shows input spectra of GPBs at various spatial frequencies ki n /k t . Gray circles 
illustrate spectrum of NRP. Lower row shows probability spectra corresponding to GPBs of input 
frequencies shown in upper row. Plots at lower part indicate azimuthal spectral distribution P(tp) at 
radius fcj n . Green and red dashed lines mark k- ln and fccBS = — k- m . 
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Figure 4. Relative power P(k) of input frequency k- m (black) and expected CBS frequency fccBS 
(blue). Spatial frequency normalized to k\ n . power normalized to total power of spectrum P tot . Errors 
are estimated by standard error of the mean power of spatial frequency. 
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